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-3 Abstract 

^' We present a QCD study of fragmentation processes for light and heavy 

^ , quarks in the semi-inclusive region of large x. Large logarithmic terms, 

due to soft-gluon radiation, are evaluated and resummed to all perturbative 
K/i ; orders in the QCD coupling as. Complete analytical results to next-to- 

;h ■ leading logarithmic accuracy are given for one-particle and two-particle inclusive 

distributions in e^e~ annihilation and DIS. Factorization of parton radiation 
at low transverse momenta is exploited to identify the universal (process- 
independent) perturbative fragmentation function that controls heavy-quark 
processes, and to perform next-to-leading logarithmic resummation of its soft- 
gluon contributions. To gauge the quantitative impact of resummation, we 
perform numerical studies of light- and heavy-quark fragmentation in e+e" 
collisions. 



CERN-TH/2001-174 
July 2001 



*This work was supported in part by the EU Fourth Framework Programme "Training and Mobihty 
of Researchers" , Network "Quantum Chromodynamics and the Deep Structure of Elementary Particles" , 
contract FMRX-CT98-0194 (DG 12 - MIHT). 

^On leave of absence from INFN, Sezione di Firenze, Florence, Italy. 



1 Introduction 



Large logarithmically-enhanced corrections due to soft-gluon radiation are a general feature 
in the perturbative-QCD study of hard-scattering processes near threshold. Techniques for 
resumming these corrections have been developed over the past several years and have been 
mainly applied to the production cross sections of high-mass systems. An extensive and 
updated list of references can be found in Sect. 5 of Ref. |]I|. 

In spite of the large amount of available data from e"^e~, lepton-hadron and hadron- 
hadron collisions (see ^^ and references therein), soft-gluon effects in single-particle (and 
double-particle) inclusive cross sections have instead received less attention. In this paper 
we consider soft-gluon resummation for the fragmentation processes of light- and heavy- 
flavoured hadrons in the vicinity of the threshold region. 

The basis for higher-order calculations in perturbative QCD is provided by the factor- 
ization theorem of mass singularities 0]. According to it, any inclusive hard cross section 
(T(a;,Q^), involving initial-state hadrons and detected final-state hadrons, can be written 
as follows 

aix,Q^)=F a D + 0{{A/QY) . (1) 

The notation in Eq. (^ is symbolic: Q^ is the hard scale, i.e. a typical transferred mo- 
mentum much larger than the QCD scale A^; x stands for any ratio of other kinematic 
invariants; the symbol ® denotes appropriate convolutions of longitudinal- and transverse- 
momentum variables and the sum over parton indices is understood. The term O {{K/QY) 
on the right-hand side of Eq. (|l|) stands for cross section contributions that are suppressed 
by some inverse power p (p > 1) of Q in the hard-scattering regime Q ^ A. 

Perturbation theory allows us to evaluate the first term on the right-hand side of Eq. (|I|), 
that is, the so-called leading-twist component of the hard-scattering process. Performing a 
power series expansion in the strong coupling as(Q^), we can compute the partonic cross 
section a and the Q^-evolution of the parton distribution functions F{x, Q"^) and of the par- 
ton fragmentation functions D{x,Q'^). Note, in particular, that only the Q^-dependence 
(or, more precisely, the anomalous dimensions) of the parton distributions -F(x, Q"^) and 
D{x,Q'^) is under control within perturbation theory. Their absolute normalization at a 
given (and arbitrary) scale has to be provided as phenomeno logical input. In the pertur- 
bative calculation this arbitrariness is reflected by the fact that anomalous dimensions and 
coefficient functions are separately dependent on the factorization scheme. Any definite 
prediction thus requires a consistent evaluation of anomalous dimensions and coefficient 
functions within the same factorization scheme. 

As long as all the kinematic scales are of the same order (i.e. x = 0{1)), perturbative 
calculations to the first few orders in the QCD coupling as should provide reliable and 
accurate theoretical predictions for the hadronic cross section. However, the perturbative 
series for the anomalous dimensions and for the coefficient functions are poorly convergent 
in the semi-inclusive or Sudakov region, that is, when the energy or transverse momentum 
of the triggered hadron (final state) is a large fraction, x — ^ 1, of the available energy ^/s in 
the scattering process. In this case the production threshold is approached and the emission 
of accompanying radiation is strongly inhibited by the kinematics. Only soft particles can 



be radiated in the inclusive final state, and the bremsstrahlung spectrum of soft (and 
collinear) gluons produces large logarithmic contributions of the type ag ln™(l — x)/{l — x) 
(with m < 2n — 1) to each order n in perturbation theory. In the presence of these 
contributions, the 'true' expansion parameter is no longer as but rather the large effective 
coupling asln^(l — x) and, hence, any finite-order perturbative calculation is unable to 
provide an accurate evaluation of the cross section. The only reliable procedure consists in 
resumming classes of logarithms to all orders in ao,. 

Leading and next-to-leading logarithmic contributions to parton distributions and par- 
ton fragmentation functions are known. To be definite, let us consider the MS factorization 
scheme and introduce the A^ moments^] Da/h,N{Q'^) of the fragmentation function of the 
parton a into the light hadron h: 



Dl%{Q')^ dxx---D^:^^M^) . (2) 



a/h,N\^ '— I ""^ -^ -^a/h 



The Q^-evolution of the parton fragmentation functions is given by the Altarelli-Parisi 
(AP) equations 



"'^a/h.N\^ _ _ 

b 



dlnQ^ 



J2^ab,NiasiQ')) Dl^^^^iQ') , (3) 



where the anomalous dimensions 7a6,Ar are the N moments of the AP probabilities. An 
important feature of the MS factorization scheme is that only the flavour-diagonal^ con- 
tributions 7qg and jgg to the evolution are affected by enhanced logarithmic corrections 
at large A^ (or, equivalently, at large x). In particular, the explicit expressions of the 
flavour- diagonal anomalous dimensions are ||5| 

7,,,^(as) ^ -Cf^ (l + K^ + 0{al))\nN +0{1) , (4) 

7,,,^(as) ^ -Ca^ {l + K^ + 0{al))\nN +0{1) , (5) 

where the coefficient K [^] is given byp| 

K=cJ'-l-^-)-'-N,, (6) 

Vi8 6 y 9 ^ ' ^ ^ 

and the term 0{1) denotes any non-singular contribution at large A^. The expressions 
(^) and d^) also show another important feature of the MS factorization scheme. In this 
scheme the anomalous dimensions are not more singular than a single power of In N when 

A^^ oo 



A similar result is valid for the parton distributions F^^^\x, Q"^), and the corresponding 
anomalous dimensions, in the large- A^ limit, are again given by Eqs. (H) and (Kf). 



•fin TV-moment space, the semi-inclusive region a; — *■ 1 corresponds to the Hmit iV — > oo. 

§The non-diagonal terms are suppressed by a relative factor of 0{1/N). 

^In this paper q;s(Q^) denotes the QCD coupling in the MS renormalization scheme. The value of the 
coefficient K reported in Eq. (g) refers to this renormalization scheme. One can also introduce an 
alternative renormalization scheme such as to absorb the coefficient K in the redefinition of as according 
to as — > as(l + Kas/'^ir). 



The knowledge of the large- iV behaviour of the anomalous dimensions is however not 
sufficient to evaluate the hadronic cross section in the large-x region. The consistency of 
the resummation procedure with the factorization theorem of mass singularities demands 
also the calculation, to the same logarithmic accuracy, of the process-dependent partonic 
cross section. The latter can strongly affect the hadronic cross section because (unlike 
the anomalous dimensions) the perturbative series for its A^ moments a^ contains douhle- 
logarithmic terms as In A^ in the general form 

n=l m=l ) 

where ajy is the leading-order (LO) contribution. Moreover, we should keep in mind that 
the soft-gluon contributions to the partonic cross section can be sizeable long before the 
threshold region in the hadronic cross section is actually approached. This is because the 
evolution of the parton densities and fragmentation functions sizeably reduces the energy 
that is available in the partonic hard-scattering subprocess; thus, the partonic cross section 
a in the factorization formula (|1]) is typically evaluated much closer to threshold than the 
hadronic cross section. 



In recent years, the general methods developed in Refs. P,P,pn}-p!^ have been applied 
to carry out soft-gluon resummation to next-to-leading logarithmic accuracy for several 
processes. Nonetheless, no explicit resummed calculation has been performed for fragmen- 
tation cross sections. 

In this paper, we consider fragmentation cross sections in e"'"e~ annihilation and deep 
inelastic lepton-nucleon scattering (DIS). We resum leading and next-to- leading soft-gluon 
contributions to one-particle and two-particle inclusive distributions. In the case of the one- 
particle distribution in e+e", we also discuss some of the next-to- next-to-leading terms. 

We also consider the fragmentation of heavy quarks. In the limit when the heavy-quark 
mass m is much smaller than the hard scale Q of the scattering process, the factorization 
formula (|I]) can be generalized in a process-independent way to compute heavy-quark cross 
sections. The generalization is based on the perturbative fragmentation function formalism 



I^,|T6|, which uses the AP evolution equations (^ to resum the large single-logarithmic 
contributions (as InQ^/m^)'^ of collinear origin. We shall show how mass and Sudakov 
effects can systematically be included in the AP evolution at low transverse momentum. 



thus extending the formalism of Refs. ||T5|,|T6[ to include soft-gluon resummation at large 



X. Our results demonstrate that in the limit m/Q -C 1, the soft-gluon contributions are 
process-independent and can thus be resummed in the universal perturbative component 
of the heavy-quark fragmentation function. In particular, we generalize the results of the 
soft-gluon resummed calculations performed in Refs. [0|l3] by evaluating the heavy-quark 
fragmentation function to next-to-leading logarithmic accuracy in the large-x region. 

The outline of the paper is as follows. In Section § we consider light-quark fragmenta- 
tion. In Sect. |2.1| we provide analytical results for the next-to-leading Sudakov resummation 



of the e^e coefficient function in the MS scheme. In the same section we also discuss some 
next-to-next-to-leading terms, and point out a universality pattern by comparing with the 



DIS structure function case. In Sect, p]^ we then perform some numerical studies and 



assess the impact of the resummation on the value of the e^e single-inclusive distribution 




Figure 1: Inclusive production of a hadron h with momentum p in e~^e annihilation. 



at large x, and on the stability of the result with respect to renormalization/factorization 
scale variations. In Section ^ we consider the case of heavy-quark fragmentation. We start 
in Sect. p.l| by briefly reviewing the perturbative fragmentation function formalism. In 



Sect. 2]2| we show how the process- independent initial condition for the perturbative frag- 
mentation function can be introduced by exploiting the universal factorization properties 
of parton radiation at low transverse momenta. Next-to-leading Sudakov resummation for 



the initial condition is performed in Sect. p.3| . The numerical studies performed in Sects. ^[3 
and ^^ assess the effects of Sudakov resummation on softening and scale dependence of the 
fragmentation function and single-inclusive e~^e~ distribution of heavy quarks. In Sect. §we 
finally summarize our main results. Soft-gluon resummation for two-particle distributions 
in e~^e~ collisions and one-hadron inclusive cross section in DIS is considered in Appendices 
A and B, respectively. 



2 Light-quark fragmentation at large x 



In this section we consider the fragmentation of light quarks (light hadrons) by performing a 
detailed theoretical and numerical study of the single-particle inclusive cross section in e~^e~ 
collisions. Related theoretical results on the two-particle distribution in e"'"e~ annihilation 
and the single-inclusive cross section in DIS are presented in the appendices. 



2.1 Single-particle inclusive distribution in e+e annihilation 

We consider the inclusive production of a single-hadron h with momentum p in e+e^ 
annihilation. Within the single- vector boson (l^ = 7*, Z^) exchange approximation (Fig. |lD, 

e+ + e- ^ViQ)^h{p)+X , (8) 

the single-particle angular distribution has three (transverse, asymmetric and longitudinal) 
components [|T^-|2D|. Each component can then be split in flavour singlet and flavour non- 
singlet contributions. In the following we do not consider the longitudinal cross section or 
the 'pure' flavour-singlet contributions, because they are suppressed by a relative factor of 



the order of (1 — x) when x ^ 1. The flavour non-singlet contributions can be written as 



da 



(e+e-), 



x,Q' 



a(^°^ /'-C(•="-)^^(x/;^,as(/i^);Q^/.^^|)D(^^)(^,/i|) , (9) 



HjJU o T 

where x is the fraction of the beam energy carried by the hadron, 

2p-Q 



^=^^' ^^°) 

and /i and y^p denote the renormalization and factorization scales, respectively. The ex- 
pression (^) is valid for both the transverse and asymmetric cross sections: a'^^'^^ are the 
corresponding cross sections at LO and D^^'^^ are the corresponding flavour non-singlet 
components of the quark fragmentation functions into the hadron h. In the following we 
consider the large-x behaviour of the (flavour non-singlet) coefficient function C'^^ ^ ^'^^(x), 
and we do not make any distinction between transverse and asymmetric coefficient functions 
because they only differ by terms that are not singular in the limit x ^ 1. 



The perturbative calculation of the coefficient functionj|j gives ||T8 



C('="-)(x,as(/.^);Q^/.^/i|) = 5(l-x) + ^^C^ 



+ 



ln(l 



X] 



X 



3 

4 VI 



X 



7r2_9 

y ~ 4 



1 /1+X2 



1 — X 

b{\ -x) + /(x 



ln^5- 



+ 0(«i) ,(11) 



where / (x) is a smooth function^ for x 
customary way: 



1 and the (. . .)+-distributions are defined in the 

dz h{z) \g{z) U^ [ dz [hiz) - hil)] g{z) . (12) 

'0 Jo 

Introducing the N moments, as in Eq. (H), the expression in Eq. (|l^) reads 

as(/i2 



C 



3 \, ,, /5 



-Cp 



IT 

9 1 



-(\nN + ^E-pi In^ + iln^iV 

4; /i^ 2 



+ (4+7.11niV+(-.^-^ + ^7l + ^7.)+0(^ 



+ 0(4 



(13) 



where 'Je = 0.5772... is the Euler constant. The large (when A^ — > oo) contributions 
In """^ N in Eq. (p!3D are due to the Mellin transformation of the singular (when x — *> 1) 
distributions [ln^(l — x)/(l — x)]+, associated to the bremsstrahlung spectrum of soft and 
coUinear emission. To higher orders in perturbation theory, the large- A^ moments of the 
coefficient function have the general double-logarithmic expansion on the right-hand side 
ofEq. (0). 

To perform the all-order resummation of the large In A^ contributions, we can apply the 
general method developed in Ref . . The single-particle distribution in Eq. (^ is obtained 

II From now on, we omit the label MS and we always use the MS factorization scheme, unless explicitly 
stated otherwise. 

**To be precise, f{x) ^ ln(l — x) when a; — > 1, and thus the contribution of f{x) to the right-hand side 
of Eq. (Ilj) is less singular than that of the other distributions in the square bracket. 



by crossing to the final state the momentum of the incoming quark in the DIS process. 
Thus, we repeat step by step and in a straightforward manner the calculation of the DIS 
cross section carried out in Ref. §]; we obtain the following resummed expression for the 
e"'"e~coefficient function: 



InC^^'''^") 



N (as(Ai );<5 ,At ,^1^) = lnAAr(as(At ),<5 //w ;(5 //Uf) 

+ In J^(as(/i'), QV/^') + Oiasias lniV)'=) 



The radiative factors An and Jn have the following exponentiated form: 



A^iasifi^), QVa*^; Q^/f^'p) = exp <^ / dz — — 



- 1 



{l-zfQ^ 



dq^ 



A[as{q')] 



IJ-F 



(14) 



(15) 



JAr(as(/i^),QVAi^ 



exp 



^N-l 



dz 



-1 



1-2 



(1~z)Q2 



dq^ 



(i-2)2q2 q 



A[as{q')] 



+ -B[as{{l-z)Q' 



n=l 



where the functions A{as) and -B(as) have perturbative expansions in as 

n=l 

whose first coefficients are: 

5(1) = -^Cf . 



- Cf ir = - Cf 
2 2 



Ca 



67_7r^ 
18 6 



Nf 



(16) 



(17) 



(19) 
(20) 



Performing the integrations over z and g^ in Eqs. (|l^) and ([T6|), we obtain a series of 
logarithmic contributions of the type ag In*" A^, with m < {n + 1). We define as leading 
logarithmic (LL) the terms with m = n+1. The next-to-leading logarithmic (NLL) contri- 
butions are those with m = n, the next-to-next-to-leading logarithmic (NNLL) terms have 
m = n — 1, and so forth. 

Owing to the knowledge of the LL coefficient A^^^ and of the NLL coefficients A^'^\ B^^^ 
in Eqs. (|19|) and (pO|), the result in Eq. (H) resums all the leading and next-to-leading IniV 
contributions to the N moments of the coefficient function of the single-particle distribution 
in e+e" annihilation. 

Note that, according to Eqs. ( [T^ ) and (|T^), in the large-A^ limit, the e^e~ coefficient 
function has the following dependence on the factorization scale fip: 



d\ndf''\ 



«s(/i^);Q^/i^/^F) 

dhifj,p 



-Alasifip)] / dz 



1 ^N-l 



(21) 



Since, to the logarithmic accuracy of Eqs. ( pjlf ) and (H), the (non-singlet) quark anomalous 
dimension 7gq,Ar(as) can be written as 

7^(as) - ^K) / dz— ^-A{as)[\nN + 0{l)] , (22) 

Jo I — z 

Eq. (^) shows that the /ii?-dependence of the coefficient function C^ ^ consistently 

matches (and thus cancels) the /ii?-dependence of the fragmentation function D)^ in 
Eq. (D. 

The NLL resummed result in Eq. (|1^ has a simple physical interpretation |]^, ^ in 
terms of independent fragmentation of the observed hadron h and of the recoiling jet. The 
radiative factor A^r describes the energy loss of the primary quark (or antiquark) that 
eventually fragments into the triggered hadron. This factor takes into account final-state 
radiation of gluons that are soft (i.e. with energy fraction uj/pq = 1 — 2;<1 — a;~ 1/A^ <C 1) 
and coUinear (i.e. with small transverse momentum g ~ a;6' = (1 — z)pq9 ^ (1 — z)Q) with 
respect to the momentum p^ of the observed hadron. Having fixed the energy of the 
jet that produces the observed hadron, the recoiling jet is constrained to carry a small 
invariant mass squared k"^ = {l — z)Q'^ ^ (1 — x)Q^ <C 1. The radiative factor J^ describes 
the fragmentation of the invariant mass of the recoiling jet as produced by the final-state 
radiation of coUinear (either soft or hard) partons. 

This independent-fragmentation picture is an effective physical picture. Although the 
two jets do not evolve independently (classically), the quantum interferences due to non- 
collinear parton radiation cancel up to NLL accuracy. The destructive interference of 
soft-parton radiation at large angles follows from QCD coherence [^, but it is no longer 
complete beyond NLL order. The resummed expression in Eq. ( p^ can be extended to any 



logarithmic order as [11-14 



Ci,^^^")(as(/i');Q',/i',/i|) = c(as(/i'),gV/^';QV^F)A!?*H«s(/i'),QV/^') (23) 

• ^Niasifi'), gV/^'; QVf^l) ^iv(as(/i'), QVf-i') + 0{l/N) . 

This expression contains two other factors in addition to those of Eq. (|14D. The factor 
c(as) does not depend on A^, and it is due to radiative corrections produced by hard (with 
energy E r^ Q) virtual partons. This factor is computable as a power series expansion in 

«s 

oo 

c(as(/i^), QV/i^; QV/il) = 1 + E ( v)" c("HQV/^^; QV/^I) • (24) 



n=l 



The radiative factor AVr° is given by 



ASr*^(as(/^'),QV/i') = exp | j^ dz ^_^ D[as{{l - zfQ')]\ , (25) 



where 



«M=(^) fl'^' + E(v)^""- P«) 

Note that A)^^ embodies In A^ contributions, and that the perturbative function -D(q;s) in 
Eq. (|26| ) has a vanishing first-order coefficient D^^\ Thus, A]^° contributes to the e"*"e~ 



coefficient function Cj^ ^ starting from NNLL order. This radiative factor takes into 
account soft-parton radiation at large angle (or, equivalently, with transverse momenta 
q r^ ijjd r^ {\ — z)Q) and it leads to violation (at NNLL accuracy) of the 'independent- 
fragmentation' picture discussed above. 

It is straightforward to check that our NLL result in Eq. ([T^ ) agrees with the large- A^ 
limit of the exact 0(«g) calculation of Ref. [^. Moreover, this calculation can also be 
used to extract a linear combination of the NNLL coefficients B^'^^ and D^'^^ in Eqs. (0) 
and (|26|). We find 

Cl {-\ + 2vr^ - 24C3) + C,C^ (-^ + ^vr^ + 40(3 



ip(2) + i5(2) ^ 1 
2 16 



+ C^^ATh'^-Itt^ 



27 9 



(27) 



where Cn is the Riemann zeta function ((^3 = 1.202 . . .). An independent calculation of -B^^-* 
and D^^-* could be performed by exploiting infrared-factorization formulae p2|p3| at 0[oi^ 
(see, e.g., the analogous calculations carried out in Refs. [0,^). The knowledge of the 
remaining NNLL coefficient A^^^ in Eq. (0) requires the (9(ag)-calculation of the anomalous 
dimensions 7gg,Ar(as) that control the evolution of the MS fragmentation functions. 

It is interesting to compare the results in Eqs. (jT^) and ( P^D with the corresponding 
resummed expression, given in Refs. 0,^, of the coefficient function for the DIS structure 
functions at large values of the Bjorken variable x. The comparison shows that, not only 
the anomalous dimensions of the parton distributions and fragmentation functions but also 
the coefficient functions exactly coincide in the semi-inclusive limit to NLL accuracy. More- 
over, the lowest-order NNLL coefficient in Eq. ( [27| ) also coincides with the corresponding 
coefficient^ for the DIS process [^. This correspondence extends, beyond the leading 



collinear level, the validity of the Gribov-Lipatov perturbative relation [^ between DIS 
structure functions and e"'"e~ fragment at ion functions. 



2.2 Numerical results 

We present some numerical results to illustrate the main quantitative effects of soft-gluon re- 
summation on the single-particle distribution in e"'"e~annihilation. Since the 0{a^ anoma- 
lous dimensions (and the NNLL coefficient A^^^') are not known, we limit ourselves to con- 
sidering NLL resummation at large-x matched to the complete next-to-leading order (NLO) 



calculation 18,19 of the coefficient function. 



Introducing the ffist two coefficients, 60 and 61, of the QCD /3- function, 

^° = VT. ' ^^ = 24^^^ ' ('^ 



^^^Note that the normahzation of our coefficients -B*^^^ and D'-^-' is different from that of the DIS coefficients 
in Ref. [^. More preciseiy, the combination (D^^^ + B^'^^)dis in Eq. (22) of Ref. |Q corresponds to our 
combination 16(D(2) + S(2)/2). 



in terms of which we haveQ 

as(/i2 



2 / A2^ 



6ilnln(^VA 



(29) 



6oln(AiVA2) V^ 6gln(/iVA2) 

and defining the variable A, 

A = 60 as (/i^) IniV , (30) 

we first evaluate the radiative factors in Eqs. (|T5|) and (p!6| ) at NLL accuracy^] and we obtain: 

lnA;v(as(/x'),QV/^';QV/^F) = lnAr/i«(A) 

+ /i(2)(A, QV/i^; QV/i|) + O (as(as IniV)^) , (31) 

lnJ^(as(/x'),QV/i') = lnAr/«(A) 

+ /(2)(A,QV/x2) + 0(as(«slnAr)'=) . (32) 

The LL and NLL functions h^^\ f^^^ and h^'^\ f^"^^ are given, in terms of the perturbative 
coefficients A^^\ A^'^\ B^^\hQM. in Eqs. (75)-(78) of Ref. 0. 

The Sudakov-resummed part Cf^ of the e^e~ coefficient function is then written (from 
Eq. (p3D) at NLL accuracy as 



+ \ --lE] In ^- 



1+^^c, 



TT 



12^ -4 + r^ + r^ 



exp 



lniV^«(A)+^(^)(A,QV/.2;QV/x|) , (33) 



where 



9 



(1) 



(A) = /^a)(A) + /«(A) 



TT^oA 



[A + (1-A)ln(l-A)] 



(34) 



^(^)(A, QV/i'; QV/i^) = h^'\\ Q'l^^ Q^l^l) + f^'-\\ QV/i') 

^ 4^ [2A + 21n(l - A) + m-d - A)] + ^^"^ -^f'^^^ ln(l - A) 



i-[A + ln(l-A)]f^-A«ln^ 



27r6o 
^Aln\. (35) 



and the term in the curly bracket is the constant (when A^ -^ 00) part of the coefficient 
function at (9(a;s) (see Eq. (|13D). 



Our final expression for the A^ moments of the coefficient function is 

cJ^^)(as(^');Q',/x',/i|) = C^(as(/x');Q',/x',/i|)- [C^(«s(/x');Q',/x',/i|)] 



as 



+ 



^{e+e ) 



c^^ ^(as(/i^);g^/i^/i^) 



(36) 



as 



where \C^ '^ j^g is the full e+e coefficient function at 0{a<s) [|l^,|T9|], C% is given in Eq. (^) 
and [C^] represents its perturbative truncation at 0(as) (i.e. at NLO). Owing to the 

*In our numerical calculations we always use the value A^^^ = 200 MeV for the QCD scale A*^^) with 
Nf — 5 effective massless flavours. This corresponds to aa,(M\^ — 0.ff6. 

^This is achieved by replacing z^^^ — f ^ —0 ( f — ^-/v~ ~ -z j in the integrand. 
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Figure 2: Effect of Sudakov resummation on a typical light-hadron fragmentation spectrum 
in e~^e~ collisions. The single-particle distribution is computed at NLO (dashed lines) and at 
NLO+NLL accuracy (solid lines). The factorization and renormalization scales are varied in the 
range Q/A < ^f = A* < 4Q. 



subtraction applied to the resummed part C^ on the right-hand side, Eq. ( ^Bf ) exactly 
reproduces the NLO result and resums soft-gluon effects beyond 0{as) to NLL accuracy. 
This defines our NLO+NLL resummed calculation. 

To obtain the single-particle inclusive cross section in Eq. (|^), the A^ moments of the 
coefficient function have to be multiplied by the A^ moments D]si{iJ?p) of the parton frag- 
mentation functions, and then we have to perfom the inverse Mellin transformation to the 
X space. Note that the resummed part of the A^ moments of the coefficient function has 
cut singularities that start at the branch-point A^ = A''^ = exp(l/6otts) (i-e. at A = 1 in 
Eqs. (^4[) and (pSf)) in the complex variable A^. These singularities, which are related to the 
divergent behaviour of the running coupling as(q'^) near the Landau pole at g = A, signal 
the onset of non-perturbative phenomena at very large values of A^ or, equivalently, when 
X is very close to its threshold value x = 1. As discussed in detail in Ref. p8|], if we are 
not interested in very high (small) values of x (Q^), we can avoid the explicit introduction 
of non-perturbative effects to deal with the Landau singularity. We thus use the Minimal 
Prescription of Ref. |2^ : the inverse Mellin transformation is performed in the complex A^ 
plane by choosing an integration contour that has all the singularities on its left, except 
for the singularity at A^ = N^, which should lie far on its right. We numerically carry out 
the inverse Mellin transformation along this contour. 

We present numerical results for the single inclusive (transverse plus longitudinal) dis- 
tribution l/cTtot da jdx in e^e~ annihilation, where the total hadronic cross section atot is 
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Figure 3: Dependence of the light-hadron fragmentation spectrum shown in Fig. ^ on the factor- 
ization/renormalization scales ^ = ^p, at three different valules of x. 



evaluated at NLO: 



0"tot(Q 



a 



(LO) 



iQ' 



1 + 



«s(g' 



TT 



+ 0{al) 



(37) 



For our illustrative purpose, we have considered only the non-singlet component of the 
distribution, and we have chosen a typical x-spectrum [|^, D{x;fj,l) = 0.11x~°'^(l — x) 
(normalized by Jq dx D{x; /ig) = 1), for the fragmentation function at the input scale /xq = 
2 GeV. The input fragmentation function is evolved up to the scale fip by using the NLO 
AP equation in the non-singlet channel (see Eq. ( ^31) for details), and then it is convoluted 
with the e~^e~ coefficient function. The (non-singlet) coefficient function is evaluated either 
at NLO or at NLO+NLL accuracy. The renormalization (yu) and factorization (fip) scales 
are fixed to be equal and are varied within the range Q/4: < fi < AQ. 

The NLO and resummed calculations at the centre-of-mass energy Q = 90 GeV are 
compared in Fig. ^j. As expected, the results of the two calculations only differ at large x. 
In the large-x region the resummed calculation has a reduced dependence on the scale /i. 
Moreover, at fixed scale fi = Q and fixed input of the fragmentation function, soft-gluon 
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Figure 4: Scaling violations at a; = 0.5 and x = 0.9, and uncertainty bands due to variations of 
the renornialization/factorization scales (/U = fip)- The fragmentation function at the input scale 
Ho is the same as in Fig. |2|. 



resummation increases the value of the single-particle distribution. This resummation effect 
can be mimicked by lowering the value of /i in the NLO calculation. 

The enhancement of the distribution in the resummed calculation can appear surprising, 
since the suppression of radiation near threshold is physically expected to decrease the cross 
section. The apparent contradiction is due to the actual definition of the fragmentation 
function and coefficient function. The effect of inhibiting soft-gluon radiation is included 
in part in the AP evolution of the fragmentation function and in part in the coefficient 
function, and the separation depends on the factorization scheme. It turns out that, in the 
MS factorization scheme, the AP evolution overestimates the physical effect of soft-gluon 
suppression. This overestimate is included in both the NLO and NLO-I-NLL calculations. 
The additional contributions that are resummed in the NLO+NLL case are those soft-gluon 
effects that are left in the coefficient function after factorization of the evolved fragmentation 
function. Having included too much soft-gluon suppression in the MS evolution, the residual 
effect in the coefficient function is positive and tends to enhance the fixed-order perturbative 
distribution. 

More numerical studies on the effect of resummation are shown in Figs. ^ and ^. 

The plots in Fig. ^ show the dependence of the fixed-order (dashed line) and resummed 
results (solid line) on the renormalization/factorization scales fi = fip, at the three points 
X = 0.1, 0.5, 0.9. As expected from Fig. |^, the NLO+NLL calculation is more stable and 
higher (when /x ~ Q) than the NLO one at large values of x. It is worth noting that the 
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agreement of the two calculations at x = 0.1 actually looks restricted to the region ^ c:i Q. 
However this statement, and in general the /i-dependence at x ^ 0.1, can strongly depend 
on the detailed shape of the input distribution D{x]iil). 

The effect of resummation as a function of the centre-of-mass energy Q is shown in Fig. ^ 
where uncertainty bands are presented for the scaling violations of the e"'"e~ fragmentation 
spectrum at x = 0.5 and x = 0.9. Lowering the energy Q, the difference between the 
resummed and fixed-order calculations is enhanced by the increase of as{Q)- The bands of 
the resummed calculation are uniformly narrower, indicating that the NLO+NLL result is 
indeed more reliable than the NLO one. 

These results on soft-gluon resummation for the MS coefficient function in e~^e~ colli- 
sions can be implemented in global fits [^ of e~^e~ data to extract parton fragmentation 
functions. In particular, the reduced scale dependence of the NLO-I-NLL calculation can 
improve the determination of the quark fragmentation functions at large x. 



3 Heavy- quark fragmentation 

In this section we consider the fragmentation of heavy quarks, i.e. quarks whose mass m 
is large enough {m ^ A) to be considered in the perturbative domain. The top quark is 
certainly heavy. Also the charm and bottom quarks are usually considered to satisfy this 
requirement. 

Heavy-quark fragmentation processes can be described by the perturbative fragmen- 
tation formalism, which resums coUinearly-enhanced perturbative contributions. In the 
large-x region, however, there are also logarithmically-enhanced contributions due to soft 
radiation. Proper resummation of both classes of logarithms is essential to ensure the re- 
liability of the theoretical predictions. In the following we first review the perturbative 
fragmentation formalism, and then discuss soft-gluon resummation. 



3.1 Collinear resummation 

The fragmentation of heavy quarks is a collinear-safe process, because perturbative collinear 
singularities are regularized by the finite value of the heavy-quark mass m. Thus heavy- 
quark fragmentation cross sections can unambiguously be computed order by order in 
QCD perturbation theory. Nonetheless, when the hard scale Q of the fragmentation 
process is much larger than m,, the perturbative series contains large logarithmic con- 
tributions, a'^{aslB-Q'^/m'^)^, of collinear origin. These logarithmic contributions have 
to be resummed to higher perturbative orders, and their resummation can be performed 
in a pro cess- independent way by using the perturbative fragmentation function formal- 
ism |lT5|,p!6|. The formalism amounts to generalizing the factorization formula in Eq. (^. 
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The cross section cfq for the inclusive production of the heavy quark Q can be written as 

^^ dz 



aQ{x,Q;m) 



/" dz 

„ Jx ^ 



(38) 
where da (o = q-, q, g) are the corresponding partonic cross sections in massless QCD and 
Da/Q is the pertubative fragmentation function of the (massless) parton a into the heavy 
quark Q. The term 0{{m/QY) on the right-hand side stands for contributions that are 
suppressed by some power p {p > 1) oi m m the kinematic regime m <^ Q. These 
contributions do not require all-order resummation and can be computed at a given fixed 
order in as- Non-perturbative corrections of the type A/Q and A/m are understood on the 
right-hand side of Eq. (|38|) . 



Note that aQ denotes a generic fragmentation cross section produced by lepton |[T5|,p!7|l 



hadron |^^|5^ or photon |^| collisions. Thus, (Tq can implicitly contain convolutions 
with parton distributions and fragmentation functions of light hadrons. The variable x in 
Eq. (|38| ) generically denotes the momentum fraction of the heavy quark {x can be either 
the energy fraction in e"'"e~annihilation or the transverse-momentum fraction in hadron 
and photon collisions). 

The resummation of the coUinear logarithms of the ratio Q'^ jvr? is achieved by writing 
the N moments of the perturbative fragmentation function as 

b 

where Eab,N{l^F, l^'op) i^ the evolution operator obtained by solving the AP equations (^, 
dEab,M,^,l,) ^ ^^^, ,,(«3(^^)) Eab,N{f^%f^lF) , (40) 

with the initial condition Eab,N{lJ^QF^ iAf) — ^ab- The starting point of the perturbative 
evolution in Eq. (|39D is set by the scale hqp, which has to be chosen of the same order as 
m, and by the initial condition E)^^}q ^, which is perturbatively computable as power series 
in as: 

oo 
^a7Q, w(«s(/io); /^0> /^0F> "^^) =5aq + Y^ as(/^o) K/Q^N(^^l^ f^OF, "^^) • (41) 

71=1 

Note that, by analogy with the partonic cross sections da in Eq. (P^, the perturbative 
expansion of the initial condition D'^Iq ^ depends on the factorization scale fiQp and on the 
renormalization scale fiQ. 



The process-dependence of the cross section (7Q(x,Q;m) is entirely embodied in the 
massless partonic cross sections da- The perturbative fragmentation function Da/Q^NifJ^'p, fn^ 
(as well as Eab,N and -D™,^ ^) is instead universal (process-independent) and can be eval- 
uated once for all. The resummation of the large coUinear logarithms is obtained by 



perturbatively solving the AP equations in Eq. ( ^OD and by setting /i^r ~ /i ~ Q and 
fiQF ~ /io ~ "^^ For instance, we can consider the second-order expansion of the flavour 
non-singlet component Jqq^N of the anomalous dimensions. 



as 



27r 



,(0) , as p(i) 



7,.,iv(as) = ;r^ Pr + 7;^Pn + 0{ai) , (42) 
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>(0) 



>(1) 



where Pjy and P^ ^^^ ^^e corresponding AP probabilities at LO and NLO. The evolution 
operator E]^{iJ?p, jj^p) for the non-singlet channel thus reads 



={0) 



ENifJ'F^ I^of) 



as(/i, 



of) 



asilA 



exp 



(as(/i^ 



of; 



"sI/^f)) 



47r26n 



P, 



(1) _ 27r6i (0) 
h ^ 



N 



(43) 






where the first factor on the right-hand side corresponds to the LO approximation, the first 
term in the curly bracket gives the NLO correction, and so forth. The leading collinear 
logarithms (as In Q"^ /m?)^ of Eq. (^) are resummed by combining the LO expression of the 
evolution operator with the LO evaluation of the partonic cross sections da and of the initial 
condition D^^Jq n- '^^^ resummation of the next-to-leading collinear terms as(cts In Q"^ /m?)^ 
requires the NLO evaluation of the evolution operator, of the partonic cross sections and 
of the initial condition. 



Although the perturbative fragmentation function is process- independent, it is not un- 
ambiguously computable. More precisely, in Eq. (|38| ) the separation between partonic cross 
sections and the perturbative fragmentation function depends on the factorization scheme. 
Since massless partonic cross sections are usually evaluated in the MS factorization scheme, 
in the following we always consider their definition in this scheme. This (implicitly) fixes 
the factorization scheme for the perturbative fragmentation function. 

Since the masses of the charm and bottom quarks are not very large, these quarks 
typically undergo a non-perturbative fragmentation process rather than perturbatively 
decay into lighter partons. In this non-perturbative process, hadronization takes place, 
eventually producing, for instance, an observable heavy meson H. Heavy-hadron frag- 
mentation cross sections can still be computed by using the factorization formula (pS]), 
provided the initial condition -D^/n for the perturbative fragmentation function is convo- 
luted with a non-perturbative fragmentation distribution D^/jj [0H131- "^^^ niost popular 
approach to describe D^m is the phenomenological model by Peterson et al. |^. Other 
approaches PB|,p7|, based on heavy-quark effective theory, are available. The hadroniza- 
tion process can also be described by modelling non-perturbative effects in the perturbative 
fragmentation function |l^,^. Recent studies |^ of non-perturbative effects may suggest 
other descriptions of D'q,^, based on a shape function to be matched to the perturbative 
part. 



At present the initial condition in Eq. (^1]) is known up to NLO. The calculation was 
first performed in Ref. [|l^] by computing the single inclusive cross section for heavy-quark 
production in e^e~ annihilation and subtracting the corresponding cross section in the 
massless case. Using the MS factorization scheme and limiting ourselves to the flavour 
non- singlet channel, the result is 



D^\x, as(/ig); ^^l dp, m') = 6{l-x) + ^^ Cp 



where the plus-distribution is defined in Eq. (|T2| 



1 + ^' fin ^0^ 
1 — X 



(1 — x)2?Ti 



2fm2 



i), 



+0(a 
(44) 
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Considering the N moments of Eq. (^) in the large- A^ hmit, we have 



TX 



2 , 2 _2n 1 , "s(AioJ 



In^ iV + ( In ^- - 27s + 1 ) In A^ (45) 






+ c?(«i) 



This expression contains logarithmic contributions proportional to In A^ and InA^. After 
having factorized and resummed large collinear terms, (as ln/x|^//iop^)", by the evolution 
of the perturbative fragmentation function, we still have to deal with the presence of In A^ 
terms in the initial condition. These terms spoil the convergence of the fixed-order per- 
turbative expansion at large A^ (or, equivalenty, at large x), and have to be resummed to 
all perturbative orders. Note that these terms are present independently of the analogous 
terms that have already been factorized at scale /xj? ~ Q and, possibly, resummed in the 
process-dependent partonic cross sections da of Eq. (|5^). The In A^ enhancement in Eq. (^) 
is a process-independent Sudakov effect due to soft-gluon radiation from the heavy quark 
at the scale /zq ~ tti. This effect can be quantitatively more important than that in the 
partonic cross sections, since it is controlled by the coupling as(Aio) ^^^^ i^ larger than 

In the following, we first discuss a general method to compute the initial condition D^^^ 
at any order in perturbation theory. We explicitly apply the method at NLO, and we 
re-derive the result in Eq. (^4]); we then perform all-order resummation of the large- In A^ 
terms up to NLL accuracy. 



3.2 Quasi-collinear factorization and the perturbative initial con- 
dition 

To understand the origin of the initial condition for the perturbative fragmentation func- 
tion, it is convenient to compare the heavy-quark cross section aQ^x, Q;rn) in Eq. (^) with 
the corresponding perturbative cross section, ag{x,Q;e), for the production of a massless 
quark q. The latter is coUinearly divergent, and thus we regularize the divergences by 
working in d = 4 — 2e space-time dimensions. In these cross sections, both the massive 
and massless quarks are produced at the large scale Q; they then perturbatively fragment 
by decreasing their momentum fraction x and radiating some amount q± of transverse mo- 
mentum. As long as qj_ is large, the massive quark behaves like a massless one and, thus, 
daQ/dq± = daq/dq^. The difference between ag and aq is produced by radiation at low 
transverse momentum, say q^ < /Iq^, where /xof is an arbitrary scale such that Q ^ yUoF 
(and /ioF > m vn the massive case). Since we are interested in the limit m/Q <^ 1, we can 
choose a value of //qf that is parametrically very small and exploit the universal factoriza- 



tion properties |40| of QCD radiation at low transverse momenta to write the A^ moments 



of the cross sections as 

(^Q,N{.Q]rn) = ^aa,NiQ;iJ'OF) D^/Q,Nil^oF,m) +0{m/Q) , (46) 

a 

fT,,^(Q;e) = ^?a,iv(Q;/^0F)5f/^,7v(/^0F,e) +0(e) ■ (47) 
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The factor a on the right-hand side conies from the high-g^ region and therefore, in the 
hmits m — ^ and e — *■ 0, it equally contributes to Eqs. (^) and (0). The fragmentation 
contributions D^m 7v(/^of, '"^) and Dj^, j^{fiQp, e), which come from the low-g_L region (q^ < 
fiQp), are different in the heavy-quark and massless- quark cases. These contributions, which 
are process-independent, are divergent in the limits m ^ and e ^ 0, respectively. In 
the heavy-quark case, the divergences are avoided by keeping m finite, although neglecting 
corrections suppressed by powers of m/fiop when m/ hqf — >■ 0. In the massless case, the 
divergences are e-poles and have to be properly factorized to obtain the short-distance 
partonic cross section aa,N that appears in Eqs. (|1|) and (pSl). The latter is defined by-^ 



0"a,7v(Q,/iOF) = ^Cr6,7v(<5;e) rjj^^^ (/iQF, c) , (4 



where ab,N{Q', e) is the dimensionally regularized cross section in Eq. (|47|) and F is a matrix 
with respect to the flavour indices a, b of the massless partons. Owing to the factorization 
theorem of mass singularities [^ , this matrix is process-independent (though factorization- 
scheme-dependent) and contains all the e-poles that cancel the singularities of a^^NiQ] e) in 
the limit e ^- 0. The explicit expression of F^'^^) as a function of the anomalous dimensions 
7afe,Ar(tts) in the MS factorization scheme is given in Ref. p. 

Inserting Eq. ( ^7|) in Eq. (^Sf ), we have a = aD^r~^. This expression for a can be 
inserted in Eq. (|38D to get erg = aD^r~^D^^\ By comparison with Eq. (^61), we thus 
obtain the final result 



7 . ^a/b,N 
b,c 



(/ioF,e) r^^^\fioF,e) D'^'Qj^{asifil);fil,filF,m'^) = D^/Q,Nif^oF,m) . (49) 

J be 



Note that the high-g_L and process-dependent contributions aa,N of Eqs. (|^) and (^) 
do not appear in Eq. (|I5|). The master equation (^) relates only process-independent 
contributions due to fragmentation at low transverse momenta, and it can be used to 
explicitly compute the initial condition D™}q j^. To this purpose we have to evaluate D^ 

and D^, cancel the e-poles in D^ by using the known expression of F*^^^), and perform the 
limits e — > and m/^Qp — > 0. 

The factorization of the low-g^ contributions D^ and D^ in Eqs. (46) and (^71) fol- 
lows from the corresponding factorization formulae (Fig. |^) of the QCD matrix elements 
(evaluated in physical gauges) . As is well known ||^, ^ ^, ^ , in the case of massless 
QCD the matrix element factorization for the splitting process a — » q{pi) + X{{qi}) {X 
denotes a set of partons with momenta qi and total momentum q = Yli qi) is controlled by 



the collinear limit. As emphasized in Ref. |^I[, the dynamics of the splitting processes of 



massive partons is described by analogous factorization formulae, which are obtained by 
generalizing the collinear limit to the quasi- collinear limit. Considering the splitting pro- 
cess a — >• Q{pi) +X{{qi}), Q being a heavy quark of mass m and momentum p'^ {pf = m^), 
we can describe its kinematics in terms of the following Sudakov parametrization: 



Pi = zp -q^ + 



z 2p ■ n 



^To simplify the notation, we set /i = ^p = Mof in Eq. (|38|). 
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Figure 5: Schematic picture of matrix element factorization in the colHnear and quasi-collinear 
limits. The thicker line denotes the particle (massless parton or heavy quark) that undergoes the 
collinear or quasi-collinear splitting process. 



z)p'' + g^ 



g^ + qi — (1 — zYrn? rf 



2p ■ n 



(50) 



denotes the collinear ('forward') direction, and 

qi = 



where the momentum p'^ (with p^ = m"^ 

the light-like vector n'^ {ri^ = 0) denotes the 'backward' direction {q±_ -p = q±-n = 0, 
—<l±q±u)- The customary Sudakov parametrization for the splitting process of a massless 
quark is recovered from Eq. (|50D by simply setting 771 = 0. We can now consider the 
collinear or quasi-collinear limits. In the massless case the collinear region is defined by 
the limit q± —* 0. In the heavy-quark case the quasi-collinear region is defined I^TJ by the 
limits q± —* and m ^ at fixed ratio m^/q5_. Note that the key difference between the 
collinear and quasi-collinear limits is that the latter has to be performed by keeping the 
mass 771 of the same order as q± while q± ^ 0. 

Performing the calculation of the collinear or quasi-collinear limits, a generic QCD 
matrix element M{pi,{qi}; . . .) factorizes in terms of a process-dependent contribution 



(which is finite in the limiting region) and a universal AP splitting function |2^, ^, ^ 
The fragmentation contributions D^ and D^ in Eq. 



19|) are obtained by integrating the 
AP splitting functions over the momenta qt of the partons involved in the corresponding 
splitting process. 

In this section we have so far outlined a general method to compute the initial condition 
jjmi £qj^ heavy-quark fragmentation. As an example, we now apply the method to the 
explicit calculation at NLO of the flavour non-singlet component of D™\ 

To this purpose, we have to consider the splitting process Q — > Q{pi) + g{q) at 0{as). 
Performing the collinear or quasi-collinear limits, the matrix element M(pi, q; . . .) factorizes 
as (see Fig. ^ 



|M(pi,g;...)p^|M(pi/^;...)|' 



Anas 
Pi -Q 



P=\Mip,/z;...)\^8nas 



zil-z) 



q^ + (1 — zYm? 



(51) 

where we have used the parametrization in Eq. ([50|), the mass-shell condition g^ = and 
the identity 

2p..,= 3l±iii:iM. ,52) 

zil - z] 
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In the massless case, P is the usual AP sphtting function in c? = 4 — 2e dimensions, 

"1 + ^2 



Pqg{z]t) 



c. 



\-z 



,{\-z) 



(53) 



while its generalization to the heavy-quark case in (i = 4 dimensions is 



PQ,(z;m7ql)=C^ 



l^z" 



m 



z V\-(l. 



C, 



l + z^ 2z{l - z)m^ 

1 — z q^ + (1 — z^m? 



(54) 



To extract the low-g^ contributions D^ and D^ in Eq. (^9]), we consider the A^ moments 
of the factorization formula (|5T|) and perform the integration of the AP splitting functions 
in Eqs. (^) and ( [5^ ) over the low transverse-momentum region, q^^ < //q^,, of the radiated 
gluon. The phase space for gluon radiation at fixed value of the longitudinal momentum 
Pi • n of the triggered quark is 



(i$ = (/i, 



\4-d 



d<^q 



^ lOTT^ 1 (1 — e) V ^1 



dz 



{2TiY-^ "^^'^ ' IGtt^ r(l-e) " "^ V qi ) z{l - z] 
where /x^ is the dimensional-regularization scale. We thus obtain 



Q{z{l-z)), (55) 



D^(/ioF,e) =1 



as 1 
27r r(l-e)jo 
as 1 1 
27r r(l-e) e 



V qi 



2\ ^ /•! 



dz iz 



N-l 



f^OF J Jo 



1) -T Piai^'^^) 
el - z) 



,(56) 



and 



1 ^"s r 



27r Jo JO 



dz iz 



N-l 



dz iz''-^ - 1 



dz [z^-^ - 1 



1 + -2^ j^ /ioF + (1 - ^)^"^^ 



q^ + (1 — zYiv? 

2z 



PQg{z-m'lcii) (57) 



A'OF 



1-2 

1 + ^2 



In- 



(1 — zYim? (1 — 2;) yUgp, + (1 — zYm? 

2z 



f^OF 



^ /yyi Zi 



zym 



+ 0{m'/iitp 



(58) 



Note that the massless- quark calculation in Eq. (^) has consistently been performed by 
using dimensional regularization, while in the heavy-quark calculation of Eq. (^) we have 
set d = A. The weight factor (2;^"^ — 1) in Eqs. ( ^6|) and ( ^8|) takes into account real 



.^iv-m 



and virtual (' — 1') gluon radiation. The real contribution derives directly from 
Eq. (^). The virtual contribution has been included by simply enforcing the constraint 
Dff^^ = -D^=i = 1, which follows from the conservation of the fermion number in the 
non-singlet sector. 



The initial condition D^j^^ is obtained from the master equation (^91), which in the non- 
singlet sector simplifies as 



^ivX"s(/io);/^0'/^OF'"^ ) 



D§{fioF,m) 



(MS), 



^^(/^0F,e) rXr 'ifioF^e) 



-1 ' 



(59) 
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or, equivalently, at 0{as): 



.(MS) 

N 



(yUoF, e) 



0(a 



(60) 

To evaluate the massless-quark contribution D^T^^, we use the C(as) expression of the 
colhnear counterterm Fj^ in the MS factorization scheme 0: 



r{MS)/ > 
Using Eqs. 



1 



"s 



1 



1 



27r r(l -e) e 
and ( 1^ we obtain 



f^OF 



2\ e /•! 





^""^'^^ ^^ziz^'-'-DCf ^ 



+ 0(ai). (61) 



-.(MS) 

N 



(/^OF, e) 



27r 



dz (z 



N-l 



1) [1 - ^] + 0(e) . (62) 



Note that this expression contains no e poles and the dependence on the dimensional- 
regularization scale fir has consistently disappeared in the limit e ^ 0. Note also that the 
finite (when e — >■ 0) term on the right-hand side entirely derives from the e-dependent part 
of the AP splitting function in Eq. (^31). 



Inserting Eq. ( |62D in Eq. (|60|) and using Eq. (|58|), we finally obtain the 0(as) contribu- 
tion to the initial condition: 



D^Ni^sif^l); f^l, I^If,^^ 



ZTT 

+ 0{al) . 



dz {z 



N-l 



l + z' 



In- 



I^OF 



zYm 



Iqm^ 



(63) 



This result agrees with the known [|T^ expression in Eq. (^4|). Our calculation explicitly 
shows that the plus-distribution in Eq. (|^) receives two contributions: the term in the 
square bracket of Eq. (^Sj) and that in the square bracket of Eq. (|62D . The former is due to 
heavy-quark fragmentation at low g^, while the latter entirely depends on the prescription 
(dimensional regularization and MS factorization) used to handle the colhnear divergences 
in the fragmentation of the massless quark. A derivation of Eq. (HI) similar to ours is 



presented in Ref. [42 



3.3 Soft-gluon resummation for the initial condition 



The 0{as) calculation presented in the second part of Sect. R]^ can be extended to the 
flavour singlet components of D^a/Q- '^^^ method discussed in the first part of Sect. ^]2| can 
also be used to evaluate D^°' at ^(ag) (and at higher orders), provided the calculation of 
the colhnear (massless) factorization formulae of Refs. [^,^ is generalized to the quasi- 
collinear (heavy-quark) case. In this section we use the master equation (|4^) (or Eq. (|59D ) 
to perform all-order resummation of the soft-gluon contributions to DJ^' at large A^ (or, 
equivalently, at large x). We limit ourselves to considering the non-singlet component, 
since the flavour- singlet contributions are suppressed by a relative factor of 0{1/N) when 
A^^ oo. 
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At high perturbative orders the fragmentation factors Dj^ and D^ receive contributions 
from muhiple radiation of partons with momenta gj. All-order resummation of these contri- 
butions is achieved as in the case of the customary AP evolution in transverse momentum. 
The region of ordered transverse momenta, fi^p > . . . > q^^j > q^^- > . . . > 0, leads to the 
exponentiation of the lowest-order AP splitting kernel. Then, the subregions where two or 
more transverse momenta are of the same order (q5_j ~ . . . ~ q^J lead to perturbative 
corrections to the exponentiated kernel. In particular, in the case of soft-gluon radiation 
the inclusive correction to the lowest-order AP kernel can be taken into account by the 
simple replacement |0,^ 



as 



as(ql) 



1 + 



as(qi) 



27r 



(K+Oie)) + Oial) 



(64) 



where the coefficient K is given in Eq. (^), and the term 0(as) only contributes beyond 
NLL accuracy at large A^. This resummation procedure applies to all the factors, Df ^ 
and D§, in Eq. (^). 



-'AT) J- N 



We first consider the case of massless-quark fragmentation. We note from Eq. (^) that 

the factor Z)^[r]y ]~^ does not show any In A^ enhancement at 0{as)- Soft-gluon radiation 
produces InA^ corrections in both Eqs. (^) and (pTj), but they cancel each other. Owing 
to Eq. (|6^) , this cancellation mechanism is valid up to NLL order in the MS factorization 

scheme, since both factors Z)^(/ioF,e) and Fjy '(fiop^e) include soft-gluon emission up to 
the transverse-momentum scale q^^ = fi^p. We thus have 



D^{^ioF,e) rf^\^oF,e) = 1 + O {asiaslnN)") 



(65) 



In the heavy-quark case, we first exponentiate Eq. (^7|) and then perform the replace- 
ment of Eq. (|64D. We obtain 



In Dj^{HoF,m) 



dz (z 



Qs(qi) 

27r 



N-l 



t^OF 



dql 



1 + ^^^K 

27r 



q^ + (1 — zYim? 
pQg{z-m'/c^l) + 0[as{as\nNf) .(66) 



The right-hand side can be further simplified by neglecting contributions of 0{m^ / ^^p) 
and terms beyond NLL order at large A^. To this purpose, we note that PQg contains 
two different contributions. In the limit m ^ 0, the q^-integrals of the first and second 
contributions in the square bracket of Eq. ( p^ are respectively logarithmically and linearly 
divergent at small q^. To the required accuracy, the logarithmically-divergent and linearly- 
divergent integrals are dominated by the transverse- momentum regions ix^p > q^^ > (1 — 
zYm? and q^ ~ (1 — zYm?, respectively. We thus obtain 



lnL'^(/ioF,H 




c?qi as(qi) 



(i-2)2m2 q± 



TT 



1 + '^^^K 



27r 



^a^{{l~zfm')\ + 0{as{as\nNf) . 



(67) 
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The expression in the curly bracket clearly shows the presence of two characteristic transverse- 
momentum scales: the factorization scale fi^p and the heavy-quark scale (1 — zYm?. The 
latter is related to the angular size 9q = m/ Eq [Eq is the energy of the heavy quark) 
of the 'dead cone' |T^ for bremsstrahlung off a massive particle. Writing the transverse 
momentum as q^ ~ (1 — z^EgO'^, where (1 — z)Eq is the energy of the radiated gluon and 
9 is its emission angle, the first term in the curly bracket of Eq. ( |67|) describes soft and 
collinear radiation outside the dead cone (^^ > 6q), while the second term is related to soft 
radiation near the dead-cone boundary {O"^ ~ ^q). 



Inserting Eq. ( pSf ) in Eq. (59), we obtain the following resummed expression for the 



initial condition of the heavy-quark fragmentation function: 

\iiD'}^\as{fily,fil,filF,m^) = \nD^{fi0F,m) + O{as{aslnN)^) , (68) 

where In D^ is given in Eq. ( |67D up to NLL accuracy. Moreover, the derivation of Eqs. (|65|) 
and ( |57D and the dependence on the factorization scale //qf of the perturbative fragmen- 
tation function in Eq. (|39|) suggest a generalization of the NLL result in Eq. ( |68D to any 
logarithmic order as follows: 

/•I r^v-i _ 1 r /■''OF Wn2 
\nD'}^\as{^^l);^^l,^^lF,m') = / dz {/ -^ A[as(qi)] 

Jo i — Z '-J(l-2)2m2 q± 

+ H[asiil-zYm^)]^ + Oil/N) . (69) 

The function A{as), whose perturbative expansion is given in Eq. (^7\), is the process- 
independent function that controls the large- iV behaviour of the flavour non-singlet anoma- 
lous dimensions (see Eq. ([2^)). The function H{as), 

oo 

i/(as)=E(v)"^'"'- '™) 

n=l 

is strictly related to soft radiation off a heavy quark, and its first-order coefficient is 

i7(i) = -Cf. (71) 

The coefficient A^^^ controls the resummation of the LL terms a^ In"^^ A^ in Eq. (|69D. The 
coefficients A^'^^ and H^^^ give the NLL terms ag In" A^, and so forth. At LL accuracy, 
Eq. (p9| ) agrees with the resummed calculation of Ref. |]T5[. Knowing the coefficients A^'^'' 



and H^^'^ in Eqs. (|^) and ([7l|), we have explicitly extended the result of Ref. |jT5[ up to 
NLL accuracy. 



To quantitatively study the effect of soft-gluon resummation, we proceed as in Sect. p72 
The Sudakov-resummed part D™' of the initial condition is written as 



D^-'iasi^,l);„l„l,,m') = {l + ^^^^ [-^ + 1 - 7I + 7i. + ( 



-^ + l-7l + 7i.+ f?-7i.Vn^ 



X exp 
where 



4 J m 

IniV^W(Ao) + (7S^(Ao,mV/ig;mV/i2^)] , (72) 

Xo = boas{fil)\nN, (73) 
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and we have explicitly introduced the renormalization scale /io, which, in general, is different 
from the factorization scale /xof- The LL and NLL functions g^J and glJ are obtained by 
performing the z and q^^ integrations of Eqs. (^) up to NLL accuracy: 



^ini(Ao) = --^[2Ao + (l-2Ao)ln(l-2Ao)], (74) 



^S(Ao,mV/i^;mV/i^^) = ^(ln^ + 27i,)ln(l-2Ao) 



[4Ao + 2 ln(l - 2Ao) + ln2(l - 2Ao 

1 /A(2) 



4nbl 



[2Ao + ln(l - 2Ao)] — - + A« In 

y-Kbo 



f4_ 

,2 



27r6o V ^^0 Kf 

+— -ln(l-2Ao). (75) 

2nbo 

Analogously to Eq. (33), the terms in the curly bracket of Eq. (^) are the constant (when 



A^ -^ oo) contributions to the NLO initial condition in Eq. (pSf ). 

The final NLO+NLL resummed expression for the N moments of the (non-singlet) 
initial condition is then given by 

T-,ini(ros)/ / 2\ 2 2 2\ 7-,ini,5/ / 2\ 2 2 2\ rnini.S/ / 2\ 2 2 2\ 

Dn ("s(/io);/^0'/^OF."^ ) = ^JV ("s(/io);^0./^OF."^ ) - [^iV ("s(/^o);/^0'/^OF."^ ) 

+ [Di^X«s(/i^);/i^,/i^^,m2)]^^ , (76) 

where [L'^'jag are the A^ moments of the C(as) result in Eq. (|^), -Dj^'' is given in Eq. (^) 
and [-D™' ]ag is its perturbative truncation at 0{as)- 

Note that, analogously to Eq. (^), the Sudakov-resummed part D]^'' of the heavy- 
quark initial condition has cut singularities in the complex variable A^. In the heavy-quark 
case the singularities start at the branch-point A^ = exp{l/26oas(A*o)} (i-^- ^^ -^o = 1/2 in 
Eqs. ( |74D and (^)). They signal the dominance of non-perturbative effects at very large 
values of A^ (or very large x). As in the massless case, we do not explicitly include any 
non-perturbative contributions in our calculations for heavy-quark fragmentation, and we 
simply apply the Minimal Prescription of Ref. ||2^ when performing numerical inversions 
of Mellin moments to x space. 

The NLO and NLO-I-NLL calculations are compared in Fig. ^ In both calculations we 
consider the perturbative fragmentation function D{x; /i|,, m"^), which is obtained from the 
corresponding initial condition D™' at the scale /xqf = Hq = m (the dependence on the 
scales Hqf and /iq is studied in Sect. p.4| ) by AP evolution (see Eqs. (^) and (^31)) up to 
the scale ^p- We can see that the NLO-I-NLL resummed fragmentation function is much 
softer than the NLO fragmentation function. Moreover, in both cases, the AP evolution 
reduces the heavy-quark momentum, and the fragmentation function is softened. 

A typical feature of all-order soft-gluon resummation is the presence of a characteristic 
maximum (Sudakov peak) [jl^ at x = Xpeak in the NLO-I-NLL fragmentation function. The 
position of the peak is at 1 — Xpeak ~ (A/mY, where the power c = 1 — exp{—7ibo/A^^^} is 
slightly smaller than unity. 
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Figure 6: Effect of Sudakov resummation for the initial condition of the heavy-quark fragmenta- 
tion function, evolved up to three different scales [Xp = 5, 50, 90 GeV. 



Figure ^ also shows that the fragmentation function becomes negative when x is very 
close to X = 1. This occurs both in the NLO+NLL calculation and in the NLO one (see also 
Fig. H), although in the NLO case the effect is not very evident from Fig. ^because it is due 
to a negative 5- function contribution concentrated near x = 1 (see Eq. (PD). In both cases 
the negative behaviour is a consequence of the approximate character of the calculations for 
the initial condition and, therefore, it is unphysical. Note, however, that the origin of this 
behaviour is different in the two cases. In the NLO case the negative behaviour is due to the 
presence of large, unresummed double-logarithmic terms as In (1 — x). In the NLO+NLL 
case this pathological behaviour is cured by soft-gluon resummation, but the sensitivity to 
non-perturbative QCD phenomena still drives the fragmentation function negative. As a 
matter of fact, if we had considered soft- photon resummation for the fragmentation of a 
massive lepton in QED, the ensuing NLO+NLL fragmentation function would have been 
positive. In QCD, instead, the Landau singularity of the perturbative coupling at small 
transverse momenta leads to a branch-point in the resummed expression ([72|) . The branch- 
point occurs when 2bQas{m'^) In A^ = 1 in the A^-plane and produces the negative behaviour 
of D{x; fijp, rn^) at 1 — x ~ K/m. At such large values of x, non-perturbative contributions 
have to properly be included in the evaluation of the heavy-quark fragmentation function 
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Soft-gluon resummation therefore suggests that non-perturbative fragmentation phe- 
nomena become dominant when 1 — x ~ K/m. This suggestion is consistent with the 
findings of the approach [§§,0 based on heavy-quark effective theory. However, the expec- 
tation of the authors of Ref. ||3^ that Sudakov effects are large only when 1 — x ~ (A/m)^, 



and likely to be less important than non-perturbative effects, is not correct. In fact, the 
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resummed calculation shows that the Sudakov effects are large at smaller values of x. In 
particular, the Sudakov peak is placed at 1 — Xpeak ~ {A/mY > A/m. Around the region 
of the Sudakov peak, perturbative soft-gluon resummation is (at least) as important as 
non-perturbative effects. 

The results presented in this section regard the process-independent fragmentation func- 
tion D{x] n'j^,m'^). This universal (though factorization-scheme-dependent) heavy-quark 
distribution enters the perturbative QCD predictions of any process-dependent cross sec- 
tions for heavy-quark fragmentation (see Eq. (pBD). In the next section, we explicitly 
consider heavy-quark fragmentation in e"'"e~ annihilation. 



3.4 Single-inclusive heavy-quark distribution in e^e annihilation 

The most accurate data on the fragmentation of bottom and charm quarks come from e~^e~ 
annihilation experiments. These data can be used to test perturbative QCD predictions and 
to extract information on the non-perturbative contribution to heavy-quark fragmentation. 

We consider the inclusive production of a single heavy meson H in e~^e~ collisions, 

e+ + e- -^V{Q)-^H{p)+X, (77) 

in the kinematical region where the centre-of-mass energy Q of the collision is much larger 
than the heavy-quark (heavy-meson) mass m. This process is completely analogous to the 
process in Eq. (^, apart from the replacement of the observed light-flavoured hadron h{p) 
by the heavy-flavoured hadron H{p). 

We study the inclusive cross section da/dx with respect to the energy fraction x = 
2p ■ Q/Q^ of the heavy-flavoured hadron, and we define the inclusive distribution 

^(x;Q^m^) = ^^, (78) 

where atot is the total hadronic cross section in e"'"e~ annihilation. Since we are mainly 
interested in the large-x behaviour, we consider only the flavour non-singlet contribution 
to the cross section. Using the perturbative fragmentation function formalism of Eqs. (^) 
and (|39|), we write the N moments of the inclusive distribution as 



5 



pAr(Q^m^) = C]v"^ ''(as(Ai^);(5^/i^/i|) SAr(/i|,/ioir) D'^\asiiJ,l); fil, filF,m 



where the ratio cr*^ '^'^V'^tot is given in Eq. (|37|) , Cj^ ^ is the MS coefficient function in 
Eq. d^), Ef4 is the AP evolution operator in Eq. (^), and D™ is the perturbative initial 
condition for the heavy-quark fragmentation function. To be precise, the right-hand side 
of Eq. ( |75[ ) should contain an additional factor D^^^^^^ describing the non-perturbative 
fragmentation of the heavy quark Q, into the observed heavy meson H . Since we are mainly 
interested in the perturbative contributions to the inclusive distribution, we do not include 
this factor. 
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Figure 7: Dependence of the moments V^ on the variations of the factorization scales ^p and 
/iOF- The renormaUzation scales /x and /io are kept fixed at Q and m respectively, except for the 
dotted lines in the right-hand plot, where /io = ^J'OF■ 



In the following we compare fixed-order and resummed calculations for the inclusive 
distribution in Eq. (^). Analogously to the comparisons in Sects. ^]2| and |3]^, the fixed- 
while the 
[TBI), which 



,{e+e-) 



©[g 



and D'^' 



order calculation uses the NLO expressions for C^ 

resummed calculation is obtained by using the expressions in Eqs. 

include the full NLO result and resum soft-gluon effects beyond 0{as) to NLL accuracy. 



and 



We recall that soft-gluon resummation for the e^e inclusive distribution V was first 

Comparing our resummed expression with the NLL results of 



considered in Ref. [jT7 



Ref. [0 in the ultrarelativistic limit Q ^ m, we find complete agreement. This should be 
regarded as a consistency check of our results. However, the authors of Ref. Jl^ did not 
use the fragmentation function formalism. Therefore, resummed formulae for the (process- 
dependent) massless coefficient function C^ '^ and the (process-independent) heavy-quark 
initial condition D}^' cannot be extracted separately from Ref. ||17|| . 

Since the most recent and accurate data on heavy-quark production in e^e~ annihilation 
come from 6-quark fragmentation at LEP and SLC, in our numerical study we choose the 
centre-of-mass energy Q = 90 GeV, and the heavy-quark mass m = 5 GeV. 

We first consider the fixed-order and resummed calculations for the N moments of the 
heavy-quark distribution ©(x; Q^, m^). Figure |^ shows the effect of varying the factoriza- 
tion scales fip and /ioF in the ranges {(5/2, 2Q} and {m/2, 2m}, respectively. The left panel 
shows the effect of varying fip, the right one that of varying fiQp. We use the customary 
practice of performing scale variations to estimate (a lower bound on) the theoretical un- 
certainty due to uncalculated perturbative terms of higher order. From the plot on the 
left-hand side, we observe that the 'theoretical uncertainty' on the large-iV moments due 
to fip variations is smaller when resummation is performed. This effect derives from the 
resummation in the e"'"e~ coefficient function Cj^ ^ of Eq. ([79|) , and thus it is very similar 



to the effect already noticed in Sect. T^ for the fragmentation spectrum of light hadrons 
(see, e.g.. Fig. |^). However, unlike those in Fig. |^, the NLO and NLO+NLL results in 
Fig. 1^ are very different, well beyond the band due to (j,f variations. The difference is due 
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Figure 8: Dependence of the moments Pat on the variation of the factorization and renormaliza- 
tion scales ^oF and /^o for a fictitious heavy quark of mass ttt, = 40 GeV. 



to the large effect produced by the resummation of the initial condition -DJ^'. Variations 
of the factorization scale /iof for the initial condition are considered in the plot on the 
right-hand side, where we can see a very remarkable reduction in the scale dependence 
of the resummed calculation. Since /iof is an auxiliary scale, introduced to perform the 
resummation of the large collinear terms as \n.Q^ /m?, the reduced scale dependence of the 
NLO+NLL calculation implies that Sudakov resummation has highly improved collinear 
resummation in the large-iV region. If we also vary the renormalization scale /xq (dotted 
lines in the right-hand plot), we observe an increased 'theoretical uncertainty' of the re- 
summed calculation. This is mainly because, in the calculation of the initial condition, we 
are using perturbation theory down to a fairly small scale /io, namely /io ~ ""^ ~ 5 GeV. 
To show how a limiting factor to the theoretical accuracy of the calculation is the low 
scale set by the bottom mass, we consider a plot analogous to the right one in Fig. |^, but 
for a fictitious heavy quark of mass m = 40 GeV. All other parameters remain identical. 
In Fig. p, we can clearly see how the 'theoretical uncertainty' is now greatly reduced by 
resummation, even when varying also the renormalization scale /xq- 

To consider also the x-space distribution, we perform a numerical inverse Mellin trans- 
formation of Eq. ([TUp by using the Minimal Prescription of Ref. [^, as already de- 
scribed in Sects. ^^ and p.3| . Figure ^ shows again the effect of varying the factorization 
scales, this time on the x-space inclusive distribution V{x;Q^,m'^). A smearing function 
jjsmear^^^ oc x°'^(l — x)^°, propcrly normalized so that L dxD^^^^^{x) = 1, has been con- 
voluted with the purely perturbative results for better visibility. As in the case of the A^ 
moments (Fig. |^, inclusion of Sudakov resummation greatly reduces the dependence on 
factorization-scale variations: the shape of the perturbative contribution to the inclusive 
distribution can now be more reliably predicted. The improvement also looks robust with 
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Figure 9: Dependence of the heavy-quark distribution ^^{x) on the variations of the factorization 
scales fip (left) and /Uqf (right). Convolution with a smearing function D^^^^^{x) oc x'^'^(l — x)'^^, 
normalized to have D^jl^^f^ = 1, has been included for better visibility. The dotted lines in the 
right-hand plot are obtained by varying also the renormalization scale /xq (/"o = I^OF = {0.5,2}). 



respect to variations of the renormalization scale /xq (dotted lines in the right-hand plot). 
We can also observe that the unphysical behaviour (its origin has been discussed at the 
end of Sect. |3]^) very close to x = 1, where the heavy-quark distribution turns negative, is 
mitigated by the inclusion of the resummed coefficient function and initial condition. 



As shown by the numerical results presented in this section (and in Sect. |0| ), Su- 
dakov resummation has important effects on heavy-quark fragmentation. In particular, 
the reduced scale-dependence of the resummed calculation permits a better control of the 
non-perturbative contributions to the heavy-quark fragmentation function. Most of these 
effects follow from the resummation in the perturbative initial condition and are, there- 
fore, process-independent. For practical phenomenological purposes, one may thus think 
of evaluating the initial condition by considering only fixed-order perturbative contribu- 
tions, without explicitly resumming Sudakov-enhanced contributions. The latter would 
then 'effectively' be reabsorbed in the definition of the non-perturbative component of the 
heavy-quark fragmentation function. However, besides it being poorly justified on theo- 
retical grounds, such an approach would also lead to much larger uncertainties due to the 
much stronger scale dependence of the fixed-order calculation. 



4 Summary 

In this paper we have performed soft-gluon resummation with NLL accuracy for fragmen- 
tation processes of light and heavy quarks (hadrons) at high momentum fraction. 

In the light-quark case, the fragmentation cross sections (see Eq. (|l|)) are obtained by 
convoluting process-dependent coefficient functions with the non-perturbative and process- 
independent fragmentation functions. We have provided the explicit resummed expressions 
for the non-singlet MS coefficient functions of the one-particle and two-particle inclusive 
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distributions in e+e" collisions and the single-particle inclusive cross section in DIS. We 
have studied in detail the e+e" one-particle distribution by matching (see Eq. (p6D) the 
NLL resummed part with the complete calculation at NLO. From the mumerical compari- 
son of our result with the NLO approximation, we observe that resummation stabilizes the 
calculation with respect to renormalization/factorization scales variations, and increases 
slightly the coefficient function in the large-x region. These findings can be useful for im- 
proved phenomenological analyses and determinations of the quark fragmentation functions 
at large x. 

In the heavy-quark case, the cross sections can be computed by convoluting coefficient 
functions, evaluated in the massless- quark approximation, with the (process-independent) 
perturbative fragmentation function of the heavy quark (see Eq. (0)). We have shown 
how the initial condition D™'(x), at the scale m, for the evolution of the perturbative 
fragmentation function can be obtained from the universal factorization properties of parton 
radiation in the collinear and quasi-collinear limits. We have then exploited this method 
to perform all-order resummation of the large soft-gluon contributions that dominate the 
behaviour of D'°'(x) at high x. We have derived an expression (see Eq. (^)) for the 
initial condition D'°'(x) that explicitly resums Sudakov terms up to NLL accuracy and 
is consistently matched with the complete NLO calculation. This result extends the LL 
process-independent calculation of Ref. []T3[ to NLL order and the NLL pro cess- dependent 
calculation of Ref. il^ in a process-independent way. 



Our numerical studies of heavy-quark fragmentation show that Sudakov resummation 
softens the high-x behaviour of heavy-quark distributions. It also leads to a marked reduc- 
tion of the dependence on the renormalization/factorization scales, making the perturbative 
predictions for heavy-quark fragmentation processes more reliable. These perturbative fea- 
tures are important for improved studies of the non-perturbative component that has to 
be included in phenomenological applications to heavy-quark fragmentation. 

Appendices 



A Two-particle inclusive distribution 
in e+e~ annihilation 

In this appendix we study the less inclusive case of e^e~ annihilation with two observed 



hadrons. Ha and Hb-, in the final state (Fig. p^Ol) : 



e+ + e- ^ V{Q) -^ hAipi) + hB{p2) + X . (80) 

We consider the corresponding two-particle inclusive distribution, 

a tot dxidx2 



1 da, . 
V{x„X2;Q')^—-^^ , (81) 
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hBiP2) 

Figure 10: Inclusive production of two hadrons with momenta pi and p2 in e~^e~ annihilation. 



where o"tot is the total hadronic cross section, and the kinematic variables Xi and X2 are 
defined^ according to Ref . |T^ : 



_ 2pi-Q _ pi ■p2 

Note that xi and X2 can independently vary in the whole kinematical range between and 
1. In the e"'"e~centre-of-mass frame, Xi coincides with the energy fraction of the hadron 
with momentum pi, while X2 is the energy fraction of the hadron with momentum p2 times 
1 — cos6'i2, 612 being the angle between the two hadrons. The choice of these two kinematic 
variables guarantees that the two-particle distribution is collinear safef] with respect to 
final-state QCD radiation for any values Xi, X2 7^ 0. Considering sufficiently large values of 
both Xi and X2, we ensure that 612 is large and, hence, the two observed hadrons belong to 
two different jets in a nearly back-to-back configuration. 

More precisely, owing to our definition of the kinematic variables Xi and X2, the two- 
particle distribution T>{xi,X2',Q'^) can be computed according to the following QCD fac- 
torization formula[] 

Vix,,X2;Q') = ^ r^ ^^C(^"^-)(a;lM,X2/^2,as(/i');Q^/i^/x|) 

<7tot Jxi ^1 Jx2 ^2 

■ D{zi,iil)D{z2^iiy) , (83) 

where D{x,^y) is the same fragmentation function as appears in Eq. (^. The two-particle 
coefficient function C^^ ^ H2;i,2;2) is computable in QCD perturbation theory. In the 
naive parton model (i.e. at the LO), we have C^'^ '^ ^(a;i,X2) = 5{1 — Xi) 5(1 — X2) and, 
thus, the two-particle distribution is simply given by the product D{xi,Q'^)D{x2,Q^) oc 
{da / dxi){da / dx2) of the two corresponding single-particle distributions in Eq. (H). Higher- 
order perturbative contributions lead to QCD correlation effects that spoil this naive fac- 
torized structure and whose size is measured by the two-particle coefficient function. 



§ Since xi and X2 are not symmetric variables with respect to pi -^^ P2, a corresponding symmetrization, 
Ha <-*■ hs, is understood in the definition of the two-particle distribution in Eq. (pl[). 

^If one uses the energy fractions 2piQ/Q^, 2p2Q/Q^ of the two hadrons in the definition of the distri- 
bution, an additional physical cutoff has to be introduced to avoid the kinematic region where pi and p2 
are collinear. 

"The ratio cr^^'-'^ /ctot is introduced on the right-hand side to make the notation in Eq. ( p3| ) consistent 
with that in Eq. (g). 
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The complete evaluation of C'^^ ^ ^(a;i,a;2) at the NLO was performed in [18|. As in 



the case of the single-particle coefficient function in Eq. (|TT]), the flavour non-singlet con- 
tributions to C^^ ^ ■'(a;i,X2) contain logarithmically-enhanced terms in the semi-inclusive 
limit xi,X2 -^ 1. These terms are conveniently singled out by introducing double Mellin 
moments as 

VnMQ^)= [ dxix^'-^ [ dx2X^'-^Vixi,X2;Q^) , (84) 

Jo Jo 



and studying their large- A^j {i = 1,2) behaviour. In A^-moment space Eq. ( p3D becomes 

riLO) 
0"tot 



VN,N,iQ') = ^:— ci;X^(as(/i');Q',/^',/iF) ^^,(/^l) D^M) , (85) 



and in the limit Aj ^ oo (i = 1,2) the coefficient function Cj^_^^^ has a perturbative 
expansion similar to Eq. (|^, with Og In"* A^ replaced by Og In"*^ A^i In™"^ A"2 (^^i +fn2 < 2n). 

All-order resummation of the logarithmically-enhanced terms can be performed by using 
standard techniques. In particular, we can exploit the strict kinematical analogy with the 
differential Drell-Yan (DY) distribution, and we can proceed as in Sect. 4 of the first paper 
in Ref. M. We obtain the exponentiated result 



fJ-F 



-1 ^A^i-l _ 1 /•(l-a.)Q^ ^^2 



Cn,n2 ("s(/i ); Q , yu , /Uf) = exp 

'0 



/ dx—^ / -\A[as{q')] + im - N2 

Jo i — X J^i q 

/ dx^ rfx2^ ^^ '-A[asiil-xi)il-X2)Q')] (86) 

Jo Jo i — Xi i — X2 

+ o(as(aslnAr,)")} . 



The structure of Eq. ( ^61) is similar to that of Eqs. (|T4D -(p!^) for the single-particle coefficient 
function. The main difference is that the integrands of Eq. (|86| ) contain a single perturbative 
function, namely the soft-gluon function A{as) given by Eq. (|l^). This is due to the fact 
that, when xi and X2 are both large, only the emission of soft gluons is kinematically 
allowed in the final state. 

The perturbative expansion of Eq. (^) at 0{as) agrees with the NLO result of Ref. [0. 
Owing to the knowledge of the first two coefficients (see Eq. (|T9D) of the function A{as), 
the result in Eq. (RBf) resums the LL (mi + 1712 = n + 1) and NLL (mi -|- m2 = n) terms 



Og In™^ A^i In"*^ A"2 in the exponent of the two-particle coefficient function Cj^^^^ . The 
integrals in Eq. ( p6D can explicitly be performed up to NLL accuracy, and the resummed 
coefficient function can be written in the following equivalent form: 



-i{e+e )/ / ,2\ ^2 ,2 ,2 



,(1)^ 



C),iA'(«s(/i');g',/,/iF) = exp[(lnAri+lnAr2) (7^;^(Ai2) 

+ ^g(Ai2,QV/^';QV/^F) + (^(«s(«slnAr,)")] , (87) 

where 

Ai2 = &oas(/i')(lnAri + lnAr2) , (88) 
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Figure 11: Inclusive production of a single hadron h' with momentum k in deep inelastic lepton- 
hadron scattering. 

and the LL and NLL functions g2p and (^^ are given by 

9i'Ji^i2) = ^^[Ai2 + (1-Ai2)ln(l-Ai2)], (89) 

^£)(Ai2,QV/i';QV/iF) = ^[2Ai2 + 21n(l-Ai2) + ln^(l-A 

- 4[Ai2 + ln(l-Ai,)]f^-A«ln^') 

- ^^V^ln(l-Ai2)-VAi2ln^. (90) 

The coefficient functions on the right-hand side of the factorization formulae in Eqs. (|^) 
and (|83|) are factorization-scheme-dependent. The resummed expressions in Eqs. ( p3D and 
(^) are valid in the MS factorization scheme. The ratio 

A(e+e-)/ / 2\ ^2 2\ ^ '-'N1N2 ('^s(/^ );Q , 1^ , I^f) /„. x 

'^k (as(/i^);Q^/i^/iF) <^iv2 (as(/i^);Q^/i^/iF) 

between two-particle and one-particle coefficient functions is instead independent of the fac- 
torization scheme. This ratio corresponds to the two-particle coefficient function as defined 
in the alternative factorization schemeQ, named 'annihilation scheme' |l^, introduced in 



Ref. ||T8[. Using our NLL resummed expressions for Cj^^ and Cj^ ^ , it is straightfor- 



ward to check that the dependence on the factorization scale hf consistently cancels in the 
ratio A}^^j^^% i.e. in the right-hand side of Eq. 



B Single-particle inclusive distribution 
in lepton— hadron collisions 



We conclude our analysis of light-hadron fragmentation in the semi-inclusive region by 
considering deep inelastic lepton-hadron scattering (DIS). Instead of studying the DIS total 

**The annihilation scheme amounts to redefining the MS fragmentation function in such a way that the 
corresponding coefficient function in Eq. (|9|) is equal to unity. 
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cross section, we are interested in the inclusive production of a single final-state hadron 
(Fig.0): 

/ + h{p) ^l' + h\k) + X . (92) 

To study the fragmentation of the tagged hadron at high momentum fraction, we can limit 
ourselves to considering the approximation in which the scattering process in Eq. (0) 
occurs through the exchange of a single vector boson with momentum q and hardness 
~(f = Q^ > 0. We denote by p and k the momenta of the incoming and final-state 
hadrons, respectively. 

Besides the customary Bjorken variable xb-, we define the final-state variable z as fol- 



lows His 



Xb = , z = , (93) 

2p ■ q p ■ q 

and the semi-inclusive region we are interested in is specified by the limit a;^ ^ 1 and 

z^l. 

In the Breit frame, where q^ = (0, 0, — Q) and p^ = {1,0,1)Q/{2xb), we have z = 
(1 — cos 6)kQ/Q, so the variable z is related to the energy feg of the fragmenting hadron 
and to its scattering angle 6 with respect to p. When approaching the semi-inclusive 
limit, the two hadrons h{p) and h'{k) are in a back-to-back configuration, with the target 
hadron momentum p'^ ~ (1, 0, l)Q/2 going forward and the fragmenting hadron momentum 
k^ ~ (1, 0, — l)Q/2 moving backward in the current-jet hemisphere. 

Note that xb and z can independently vary in the whole kinematical range between 
and 1. Moreover, as long as z is not vanishing, the fragmenting hadron h'{k) cannot 
become coUinear to the colliding hadron h{p). Therefore, the single- inclusive cross section 
daf^, /dxBdz fulfils |]TB|] the following QCD factorization formula: 

d'^h' (0) [^ dx f^ dy 



aW 



Jxr, X J, y 



dxB dz 

■ F{xB/x,fj:'p)Diz/y,ij:j,) , (94) 



.2 \ n/-,/„. ,.2 



where a^^' is the LO cross section, D(z,jj,'jp) is the parton fragmentation function of the 
tagged hadron h'{k) and F{x, /i|,) is the parton distribution function of the colliding hadron 
h{p). To be defined, we recall [Q that F(a;, /i|,) appears in the analogous factorization 
formula for the DIS total cross section da^^^^^ /dxB'- 



= cr^ 

dXB 



r dr 

/ — c^^'^\x,asi^i')■,Q^^i^^^l) FixB/x,fil) . (95) 

JxB -^ 



Since eventually we are mainly interested in the semi- inclusive limit , we omit parton indices 
in Eqs. ( |9^ ) and (pSj), and we understand that the parton fragmentation function and 
distribution function refer to the flavour non-singlet components. 

The single-inclusive coefficient function C^^^^''{x,y) is computable in QCD perturbation 
theory. In the naive parton model (i.e. at the LO), we have C-^^^-* (x, y) = 5(1 — x) 5{1 — y), 
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and the single-inclusive cross section in Eq. ( pi]) is simply proportional to the product of the 
parton distribution and the fragmentation function, daj^, /dxBdz oc F{xb,Q'^)D{z,Q'^). 
The complete NLO calculation of C^^^^\x, z) was performed in Ref. |T^ 



As in the case of the e~^e~ two-particle coefficient function C^^ ^ •'(a;i,a;2) in Eq. 
at high perturbative orders the flavour non-singlet contribution to C^^^^\xb, z) contains 
terms that are logarithmically enhanced in the semi- inclusive limit xb,z -^ 1. Similarly to 
the two-particle distribution in e~^e~ annihilation, it is thus convenient to define the double 
Mellin moments: 

cS(«s(/^');Q',/^',/i|) = Trfxx^-i r dz z''--'C^'''^\x,z,as{fi'y,Q',fi',fil) , (96) 

Jo Jo 

and to consider their large- A^j {i = 1, 2) behaviour. 

The general structure of C'jy^jVa i^ ^^^ large- Aj limit is similar to that of C'jy^^^ ^^^ 
the resummation of the large logarithmic contributions can be performed following Ref. ||^. 
Using the MS factorization scheme, to NLL order we find 

lnC'S(«s(/w');Q^/i^/i|)=lnC'J,^X^(as(/i2);Q^/x^/i?.) + 0(«s(^^ , (97) 

with Cj^_^^^ given by the expression in Eq. ([861) . This result has a simple physical explana- 
tion. In the semi-inclusive limit, the kinematical configuration of the DIS single- inclusive 
cross section is related to that of the e~^e~ two-particle cross section by crossing the DIS 
incoming hadron to the final state. In both cases the cross section is dominated by purely 
soft emission. The relation ( p7| ) thus follows from the fact that the intensities of soft-gluon 
radiation from space-like and time-like partons are equal up to NLL accuracy (see Ref. [0 
and the discussion at the end of Sect. |2?T| ). 
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